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1.  Overview 

Studies  of  elementary  reactions  that  are  of  importance  in  chemical  iodine  lasers  were 
the  primary  focus  of  the  work  described  in  this  report.  One  motivation  for  these  studies 
was  provided  by  the  need  to  achieve  a  better  understanding  of  the  kinetics  of  Chemical 
Oxygen  Iodine  Laser  (COIL).  The  pump  reaction  of  this  device, 

02(a1A)  +  I(2P3/2)  02(X35T)  +  I(2P1/2)  (1) 

was  studied  at  low  temperature  (T=150  K)  to  provide  rate  constant  data  suitable  for 
modeling  supersonic  lasers1’2.  This  reaction  was  also  investigated  using  high  level 
theoretical  methods,  in  order  to  test  the  ability  of  such  models  to  predict  electronic  energy 
transfer  kinetics3  4. 

Iodine  lasers  operate  on  a  single  hyperfine  of  the  transition 

(F’=3— >F”=4).  As  a  consequence,  power  extraction  is  influenced  by  the  rate  at  which 
collisions  transfer  population  between  the  F’=2  and  F’=3  sub-levels  of  2P m.  The 
hyperfine  transfer  rate  constants  had  not  been  measured  previously,  and  most  models  of 
COIL  systems  have  assumed  that  the  process  is  fast  enough  to  maintain  equilibrium 
populations  in  the  F’=2  and  3  sub-levels.  In  the  present  work  we  performed  the  first  direct 
measurements  of  the  hyperfine  transfer  rate  constants  over  the  temperature  range  from  10  - 
300  K1’2.  The  rate  constants  were  smaller  than  had  been  assumed,  indicating  that 
computational  simulations  of  COIL  should  include  specific  treatment  of  hyperfine 
relaxation. 

Characterizations  of  a  prototype  supersonic  COIL  device  (VertiCOIL)  revealed 
cavity  temperatures  that  were  appreciably  greater  than  the  value  predicted  by  computational 
fluid  dynamics  models5.  It  was  suggested  that  the  heat  release  was  due  to  either 
condensation  of  HjO  in  expansion  or  quenching  of  I(2P1/2)  by  F^O.  Experiments  designed 
to  test  both  of  these  possibilities  were  carried  out.  The  results  did  not  support  either 
hypothesis. 

The  mechanism  by  which  I2  is  dissociated  by  02(a‘A)  is  one  of  the  least  well 
understood  aspects  of  COIL  systems6.  It  has  been  estimated  that  somewhere  between  four 
and  six  02(a)  molecules  are  needed  to  dissociate  one  I2  molecule,  but  the  details  of  the 
process  remain  obscure.  We  have  investigated  kinetic  models  of  the  dissociation  process 
that  identify  vibrationally  excited  I2(X)  as  the  key  intermediate  species7.  We  were  unable 
construct  satisfactory  models  based  on  this  mechanism,  and  have  proposed  a  modified 
scheme  that  involves  electronically  excited  I2  (a  variant  of  models  originally  proposed  by 
Bacis  and  co-workers8).  Measurements  of  quenching  rate  constants  for  the  A’  state  of  I2 
were  carried  out  to  obtain  parameters  for  this  model.  The  revised  mechanism  provides  a 
better  quantitative  description  of  the  dissociation  kinetics. 
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Iodine  lasers  that  are  driven  by  energy  transfer  from  NCl(a)  are  currently  being 
developed911.  Measurements  of  the  rate  constants  for  quenching  and  reaction  of  NCl(a) 
are  needed  to  guide  the  design  and  development  of  such  devices.  Earlier  studies  of  NCl(a) 
quenching  yielded  rate  constants  that  were  strongly  dependent  on  the  technique  used  to 
perform  the  measurement1214.  In  some  instances  different  techniques  yielded  rate  constants 
that  differed  by  more  than  an  order  of  magnitude.  New  measurements  were  performed  to 
resolve  these  conflicts  and  provide  data  for  previously  unexamined  collision  partners15. 
Self-annihilation  of  NCl(a)  was  one  of  the  most  significant  processes  examined.  This 
second  order  reaction  limits  the  concentration  of  NCl(a)  that  can  be  generated  by  chemical 
means  and  dominates  the  transport  losses.  We  obtained  a  self-annihilation  rate  constant  of 
(7±2)xl0‘13  cm3  s'\  which  is  a  factor  often  smaller  than  the  previous  estimate13.  This  result 
is  encouraging  for  NCl(a)  transfer  lasers  as  it  indicates  that  transport  and  mixing  losses  will 
be  less  sever  than  earlier  kinetic  models  had  predicted. 

Quenching  of  NCl(a)  by  I2  was  found  to  be  a  reactive  process  that  liberates  atomic 
I.  Subsequently,  the  electronic  energy  transfer  process  I(2P3/2)+NCl(a)  — >  I(2P1/2)+NCl(X) 
lead  to  emission  from  atomic  iodine.  Analysis  of  the  temporal  behavior  of  the  iodine 
fluorescence  indicated  that  quenching  of  I(2P1/2)  by  C1N3  (the  NCl(a)  precursor)  was  very 
slow.  This  was  confirmed  in  a  study  of  the  quenching  of  I(2P1/2)  by  C1N3  and  HN3. 

A  method  for  observing  NC1(X)  by  means  of  the  b-X  absorption  bands  was 
developed15.  This  is  of  interest  as  absorption  measurements  may  be  used  to  characterize  the 
temperature  and  pressure  conditions  within  the  cavity  of  an  NCl(a)/I  laser  system.  The 
absorption  measurements  were  carried  out  in  collaboration  with  Dr.  G.  C.  Manke  II  (AFRL, 
NM)  and  Dr.  S.  J.  Davis  (PSI,  MA).  The  results  are  being  used  in  the  design  of  a 
diagnostic  instrument  for  the  AFRL  azide  laser  flow  facility.  Time-resolved  NC1(X) 
absorption  measurements  were  used  to  determine  the  yield  of  NCl(a)  from  C1N3  photolysis. 

Theoretical  studies  of  I+NCl(a)  energy  transfer  and  the  reaction  Cl+N3 
— »NCl(a,X)+N2  were  initiated  in  a  collaborative  effort  with  Professor  K.  Morokuma.  The 
main  objective  of  these  projects  is  to  determine  product  branching  ratios.  The  efficiency  of 
an  NCl(a)/I  laser  system  is  dependent  on  these  branching  ratios,  which  are  difficult  to 
measure  accurately  by  experimental  means.  For  I+NCl(a)  the  potential  energy  surfaces 
have  been  calculated,  and  the  surface  intersections  responsible  for  electronic  energy  transfer 
have  been  located.  As  the  relevant  potentials  do  not  correspond  to  chemically  bound  states, 
it  appears  that  energy  transfer  occurs  by  a  purely  physical  mechanism.  Calculations  for  the 
Cl+N3  reaction  have  been  used  to  locate  the  barriers  to  formation  of  singlet  (NCl(a)+N2) 
and  triplet  (NC1(X)+N2)  products.  Both  barriers  lie  below  the  energy  of  the  entrance 
channel.  From  the  energetics  alone  we  cannot  explain  the  observed  preference  for 
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formation  of  the  singlet  product.  Dynamical  calculations  are  planned  that  will  probe  this 
question  in  more  detail. 

2.  Studies  of  I(2Piy2)+02  Energy  Transfer. 

2a.  Experimental  studies  of  electronic  energy  transfer  and  nuclear  spin  state  changing 
collisions. 

The  energy  transfer  process 

02(a‘A)  +  I(2P3/2)  02(X32T)  +  I(2P1/2)  (1) 

is  exothermic  by  0.8  kcal  mol'1  (E/R=401.4  K)  in  the  forward  direction.  For  the  purpose  of 
modeling  COIL  kinetics,  rate  constants  for  both  the  forward  and  backward  transfer  steps  are 
needed.  In  practice,  the  rate  constant  for  the  reverse  reaction  (kb)  is  most  easily  measured. 
It  is  known  that  the  quenching  of  I*  (=  I(2P1/2))  by  02(X)  is  dominated  by  the  E-E  energy 
transfer  process.  Hence,  measurements  of  the  quenching  rate  constant  yield  good  estimates 
for  k,,.  Once  k,,  is  known,  the  forward  rate  constant  (lq )  can  be  calculated  reliably  from  the 
theoretical  expression  for  the  equilibrium  constant.  Rate  constants  for  the  forward  and 
backward  transfer  at  room  temperature  are  well  established.  However,  for  the  purpose  of 
modeling  supersonic  COIL  devices,  knowledge  of  the  temperature  dependence  of  the  rate 
constants  over  the  range  150<T<300  K  is  needed. 

In  this  program  we  characterized  I*  +  02  energy  transfer  at  temperatures  near 
150  K.  As  in  our  previous  work  on  this  system16,  a  Laval  nozzle  expansion  was  used  to 
generate  low  temperature  gas  flows.  I*  was  produced  by  pulsed  photolysis  of  I2.  Two 
different  time-resolved  techniques  were  used  to  determine  the  rate  constant  for  energy 
transfer  from  I*  to  02.  Initially  a  fast  IR  detector  was  used  to  observe  I*  fluorescence  decay 
curves.  Measurements  of  the  decay  rate  as  a  function  of  02  concentration  yielded  a  transfer 
rate  constant  of  kb(150)=(6.3±1.0)xl0’12  cm3  s’1.  This  value  was  somewhat  higher  than  our 
previous  determination16  of  kb(150)=(4.6±1.0)xl0’12  cm3  s’1  (measurements  that  involved 
deconvolution  of  signals  from  a  slow  Ge  detector). 

Pulsed  LDF  detection  of  I*  was  used  for  the  second  series  of  measurements.  Our 
initial  attempts  with  this  method  yielded  an  energy  transfer  rate  constant  of 
kb(150)=(7.6±1.0)xl0'12  cm3  s’1.  We  speculated  that  the  discrepancy  between  the  rate 
constants  derived  from  LIF  and  IR  fluorescence  may  be  indicative  of  systematic  errors  in 
the  former.  Further  investigation  of  this  discrepancy  revealed  a  rather  surprising  aspect  of 
I2(J3)  photolysis. 
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In  the  LIF  experiments,  I*  was  detected  by  exciting  the  5p46s  (2P3/2)  -  5 p5  (2P m) 
transition.  Owing  to  the  large  hyperfine  splitting  of  the  2Pm  state,  the  spectrum  shows  two 
resolved  features,  corresponding  to  transitions  from  the  I*  F=2  and  F=3  sub-levels.  Fig.  1 
shows  spectra  taken  using  three  different  methods  for  production  of  I*.  The  upper  and 
lower  traces  correspond  to  photolysis  of  I2  at  498  nm  and  490  nm,  respectively.  The  center 
trace  was  obtained  when  I*  was  generated  by  248  nm  photolysis  of  CF3I.  It  is  apparent 
from  these  traces  that  the  different  photolysis  sources  yielded  different  initial  populations  in 
the  atomic  hyperfine  sub  levels.  Photolysis  of  CF3I  is  known  to  result  in  statistical 
population  of  the  I*  hyperfine  levels  (i.e.,  [F=3]/[F=2]  =  7/5).  This  distribution  is  evident 
in  the  center  trace  of  Fig.  1.  The  upper  trace  shows  that  498  nm  photolysis  of  I2  results  in 
preferential  population  of  F=2. 

When  quenching  of  I*  by  02  was  investigated  using  498  nm  photolysis,  the  kinetics 
of  the  F=2  and  F=3  sub  levels  appeared  to  be  different.  For  example,  Fig.  2  shows  the 
effect  of  adding  02  to  the  expansion.  This  trace  gives  the  impression  that  F=2  is  quenched, 
while  F=3  is  not  affected.  It  is  most  unlikely  that  the  hyperfine  sub-levels  are  quenched 
with  different  rate  constants.  The  most  probable  explanation  for  this  behavior  is  that  the 
quenching  rate  constants  for  the  hyperfine  sub-levels  are  the  same,  but  the  kinetics  are 
complicated  by  transfer  between  the  hyperfine  levels  (initially,  transfer  from  F=2  to  F=3). 
Kinetic  modeling  of  the  data  supported  this  interpretation.  Analysis  of  the  LIF  data  yielded 
an  electronic  energy  transfer  rate  constant  of  kb(150)=7.3xl012,  and  hyperfine  transfer  rate 
constants  of  k(2->3)=2.5xl0n  and  k(3-»2)=1.8xl0u  cm3  s’*.  In  accordance  with 
theoretical  predictions,  we  found  that  the  closed-shell  collision  partners  He,  Ar,  and  N2  did 
not  induce  measurable  degrees  of  hyperfine  relaxation  under  the  conditions  of  our 
experiments. 

For  I*+02(X)  an  optimized  fit  to  both  the  IR  fluorescence  and  LIF  data  gave  rate 
constants  of  1^(1 50)=(7.0±0.7)xl012,  k(2-43)=2.6xlO  n  and  k(3->2)=1.9xl0'H  cm3  s'1. 
The  revised  value  for  1^(150)  was  significantly  below  the  rate  constants  predicted  by 
currently  accepted  models.  When  combined  with  the  room  temperature  value  for  1^,  the 
present  result  is  consistent  with  the  Arrhenius  expression 

kb(T)  =  1.05xl0"‘°  exp(-401.4/T)  cm3  s1 

The  rate  constant  for  02(a)+I  ->  02(X)  +  I*  was  derived  from  this  expression  via  the 
statistical  mechanics  expression  for  the  equilibrium  constant.  A  temperature  independent 
value  ofkl=7.8xl011  cm3  s”‘  was  obtained.  Further  studies  were  carried  out  to  characterize 
the  temperature  dependence  of  the  I*(F)+02(X)  — »I*(F’)+02(X)  hyperfine  transfer 
process.  To  augment  energy  transfer  data  taken  at  150  K,  measurements  were  made  at 
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10  K,  in  a  free-jet  expansion,  and  295  K  in  a  conventional  vacuum  system.  As  part  of  this 
investigation  we  characterized  the  nascent  I*  hyperfine  level  population  distributions 
resulting  from  photolysis  of  Ij  at  different  temperatures  and  wavelengths. 

Nascent  I*  population  distributions  resulting  from  photolysis  of  jet-cooled  I2  (10 
K)  were  examined  for  several  photolysis  wavelengths  in  the  490-499  nm  range.  Examples 
of  the  I*  I  IF  spectra  are  shown  in  Fig.  3.  As  the  photolysis  wavelength  increased  from  490 
to  499  nm,  the  probability  for  populating  the  F=3  level  decreased  dramatically.  Note  that 
499  nm  corresponds  to  the  onset  of  the  B-X  continuum.  The  lower  trace  in  Fig.  3  shows  that 
dissociation  at  energies  just  above  the  B  state  dissociation  threshold  strongly  favored 
formation  of  F=2  products.  Population  transfer  between  hyperfine  levels  is  most  easily 
observed  when  the  initial  population  distribution  is  far  from  equilibrium.  Consequently,  499 
nm  photolysis  was  used  for  energy  transfer  measurements.  Fig.  4  shows  spectra  taken  in 
the  free-jet,  with  and  without  02  present.  Comparing  these  traces,  it  can  be  seen  that 
collisions  with  02  are,  on  average,  moving  population  from  F=2  to  F=3.  At  the  local 
temperature  of  these  measurements,  the  rate  constant  for  electronic  quenching  of  I*  by  02 
should  be  negligibly  slow  (the  main  channel,  transfer  to  02(a),  is  exothermic  by  278  cm  *). 
Analysis  of  these  data  yielded  rate  constants  of  k(2— >3)=2.0±0.5xl0 11  and 
k(3-*2)=1.4±0.5xl0"  cmV  at  T=10  K. 

Nascent  I*  distributions  resulting  from  \  photolysis  at  room  temperature  woe  also 
non-statistical,  but  the  effect  was  not  as  pronounced  as  it  was  for  photolysis  at  10  or  150  K. 
As  for  the  low-temperature  experiments,  the  relative  population  of  F=3  decreased  with 
increasing  wavelength.  At  the  long  wavelength  limit  (499  nm)  the  [F=3]/[F=2]  population 
ratio  was  approximately  2/3.  Although  electronic  quenching  was  as  fast  as  hyperfine 
transfer  at  room  temperature,  the  latter  was  still  readily  observable.  Kinetic  modeling  of  the 
room  temperature  data  defined  rate  constants  of  k(2— »3)=(2.3±0.5)xl0 11  and 
k(3-»2)=(1.6±0.5xl0'“)  cm3 s'1.  The  hyperfine  transfer  rate  constants  were  found  to  have 
no  discemable  temperature  dependence  over  the  range  10-295  K.  As  the  hard-sphere 
collision  frequency  decreases  with  temperature,  this  suggests  that  the  transfer  cross  sections 
increase  with  decreasing  temperature. 

For  computational  models  of  COIL  systems  it  is  usually  assumed  that  hyperfine 
transfer  is  a  rapid  process.  The  relative  populations  in  the  F=2  and  F=3  sub-levels  are  fixed 
at  the  equilibrium  value,  and  it  is  assumed  that  the  energy  stored  by  I*(F=2)  is  available 
through  hyperfine  transfer.  Our  measurements  show  that  the  hyperfine  transfer  rate 
constants  are  smaller  than  theoretical  estimates,  and  that  the  equilibrium  assumption  is  not 
justified. 
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Figure  3.  I*  LIF  spectra  for  photolysis  of  I2  at  490, 498,  and  499  nm. 
These  data  were  taken  at  10  K  in  an  He  expansion  with  a 
photolysis  -  probe  delay  of400  ns. 
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Figure  4.  I*  LIF  spectra  for  photolysis  of  I2  at  499  nm  in  a  ffee-jet  expansion. 
The  upper  trace  shows  the  nascent  population  distribution  of  the  hyperfine 
sub-levels.  The  lower  trace  shows  the  effect  of  adding  02  to  the  expansion. 
Both  traces  were  recorded  using  a  photolysis-probe  delay  of  800  ns. 


10 


2b.  Comment  on  the  photodissociation  dynamics  ofI2(B). 

The  non-statistical  population  of  I*  hyperfine  levels  produced  by  IfB)  photolysis 
was  an  unexpected  result.  In  the  absence  of  electronic  state  perturbations,  dissociation  of 
I2(B)  should  yield  a  statistical  distribution.  TJie  observed  distributions  were  consistent  with 
electronic  state  mixings  that  are  mediated  by  the  J-f  term  in  the  molecular  Hamiltonian. 
This  operator  can  mix  states  that  differ  in  Q.  by  ±1,  and,  in  the  process,  break  the  u/g 
symmetry  of  the  molecule.  Such  symmetry-breaking  perturbations  near  the  1 2(B) 
dissociation  limit  have  been  observed  previously17.  Avoided  crossings  between  curves  that 
correlate  with  I*(F=2)+I  and  I*(F=3)+I  can  influence  the  final  state  population 
distribution,  provided  that  the  atoms  move  relatively  slowly  through  crossing  regions 
(adiabatic  curve  following).  The  fact  that  the  lower  energy  hyperfine  level  receives  the  most 
population  is  consistent  with  an  avoided  crossing  model.  Note  that  the  strong  preference 
for  population  of  F=2  seen  for  499  nm  photolysis  of  the  10  K  sample  (lower  trace  in  Fig. 
3)  is  not  simply  an  energetic  threshold  effect.  The  photon  energy  is  above  the  I*(F=3)+I 
dissociation  limit. 

As  the  velocity  at  which  the  I  atoms  traverse  the  curve  crossings  increases,  the 
dynamics  switch  from  adiabatic  to  diabatic  curve  following.  Hence,  as  the  photolysis 
wavelength  is  decreased,  the  hyperfine  population  distribution  moves  toward  the  statistical 
limit.  This  trend  is  evident  in  Fig.  3.  For  a  fixed  photolysis  wavelength,  the  nascent 
distribution  clearly  exhibits  a  dependence  on  temperature.  As  the  temperature  increases,  the 
[2]</[3]0  ratio  decreases.  This  behavior  may  also  be  explained  by  a  curve  crossing  model. 
Thermal  energy  simply  adds  to  the  energy  imparted  by  the  photon,  and  thereby  increases 
the  final  recoil  velocity.  Consequently,  the  [2]</[3]0  ratio  observed  for  495  nm  photolysis  at 
10  K  was  very  similar  to  that  seen  for  498  nm  photolysis  at  150  K. 

In  future  experiments  it  will  be  of  interest  to  probe  the  nascent  distribution  of  1 2P3/2 
hyperfine  levels  populated  by  I2  B  state  photolysis.  The  distributions  among  these  levels 
may  also  be  non-statistical.  If  this  proves  to  be  the  case,  near-threshold  photolysis  of  I2(fl) 
may  provide  a  convenient  method  for  studying  hyperfine  relaxation  of  ground  state  atoms. 

2c.  Theoretical  calculation  of  the  rate  constants  for  I(2P3/2)  +  02(a'Ag)  <=>  I(2PI/2)  + 
02(X3£g)  ( collaboration  with  Prof.  K.  Morokuma). 

Theoretical  studies  of  the  l+02(a)  energy  transfer  process  were  performed  to  obtain 
a  deeper  understanding  of  the  factors  that  determine  the  temperature  dependence  of  the  rate 
constant.  This  system  also  provides  a  challenging  test  of  current  theoretical  methods. 
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Two  dimensional  potential  energy  surfaces  for  I(2P1/2)  +  02(X)  and  I(2PM)  +  02(a) 
were  calculated  using  high-level  ab  initio  methods  (CASSCF  followed  by  CASPT2  with  the 
cc-aug-pvtz  basis  set)4.  The  0-0  bond  was  held  fixed  at  the  centroid  value.  The  Jacobi 
coordinates  that  describe  the  position  of  the  I  atom  relative  to  the  0-0  bond  (R  and  y)  woe 
varied  to  generate  a  grid  of  energy  points  for  each  electronic  state.  Seven  potential  surfaces 
were  obtained  and  then  fitted  exactly  with  simple  analytical  functions  of  the  type 

v:j(r.R)=  'icf(R)p,(cos(r))  (2) 

1=0,2,... 


where  cf  (R)  are  radially  dependent  expansion  coefficients  interpolated  using  1-D  splines, 


and/^cos^))  are  normalized  Legendre  polynomials.  The  indices  i  and  j  denote  the 

electronic  state.  Diabatic  coupling  elements  connecting  different  pairs  of  states  i^j  were 
calculated  on  the  same  grid,  and  fitted  with  the  same  type  of  analytical  functions3. 

To  a  good  approximation,  we  can  treat  the  nuclear  dynamics  classically  since  we  are 
dealing  with  interactions  of  heavy  particles.  Hamilton’s  equations  of  motion  were  solved  to 
propagate  the  position  of  the  nuclei  in  time  according  to  the  momenta,  P  (radial)  and  j  (02 


rotation). 


p2  +-^+tLJt+Vu(r,R) 


2ju  2  fir*  2  nR 


(3) 


Eq.  (3)  is  the  Hamiltonian  for  the  k-th  potential,  where  jU  is  the  reduced  mass  of  the 


collision  complex,  fX  is  reduced  mass  of  02  and*J  is  the  total  rotational  angular 
momentum.  The  corresponding  equations  of  motion  for  coordinates  and  conjugate 
momenta  were  integrated  using  a  fourth  order  Runge-Kutta  algorithm. 

The  electronic  degrees  of  freedom  must  be  treated  quantum  mechanically.  For  these 
calculations  we  used  the  surface  hopping  trajectory  methods  developed  by  Tully18.  Electron 
densities  were  propagated  via  the  time-dependent  Schrodinger  equation: 


daki 

ih — - 
dt 


(4) 


where  the  complex  coefficients  akj  are  amplitudes  whose  squared  diagonal  elements  are 

populations,  and  off-diagonal  elements  are  coherences.  The  nuclei  and  electron  densities 
were  propagated  simultaneously,  and  transition  probabilities  were  calculated  according  to 
Tully’s18  “fewest  switches”  model. 
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To  calculate  the  temperature  dependence  of  the  rate  constant  for  I(2P3/2)  +  02(a‘Ag) 

I(2P1/2)  +  02(X3Zg')  approximately  5x10s  trajectories  were  statistically  averaged  over  a 
range  of  impact  parameters  and  collision  energies.  Temperatures  in  the  range  10-300  K 
were  investigated.  It  was  found  that  collisions  resulting  in  energy  transfer  were  dominated 
by  single  hop  trajectories.  The  calculated  energy  transfer  rate  constant  decreased  smoothly 
with  increasing  temperature  over  the  range  100-300  K.  The  predicted  value  was  in  excellent 
agreement  with  the  experimental  result  for  150  K,  but  the  calculations  underestimate  room 
temperature  data  by  a  factor  of  1.6.  The  rate  constant  increases  with  decreasing  energy 
because  (i)  long-range  attractive  forces  draw  slow  moving  collision  partners  together  and, 
(ii)  the  longer  duration  of  slow  collisions  increases  the  probability  of  surface-hopping.  It 
was  also  found  that  there  is  a  competition  between  rotational  relaxation  of  O 2(a)  and 
electronic  energy  transfer. 

Given  the  complexity  of  this  system  we  considered  that  the  level  of  agreement 
between  theory  and  experiment  was  very  encouraging.  The  causes  of  the  residual 
discrepancies  could  be  related  to  errors  in  the  potential  energy  surfaces  and  neglect  of 
quantum  interference  effects. 

4.  Measurements  related  to  heat  release  due  to  H20  in  COIL  systems. 

Studies  of  the  VertiCOIL  laser  system  at  Kirtland  AFB  revealed  unexpectedly  high 
cavity  temperatures5.  Model  calculations  predicted  downstream  temperatures  of  160  K, 
whereas  the  measured  values  were  around  220  K.  This  was  a  matter  of  concern  for  the 
groups  who  simulate  COIL  systems,  as  the  discrepancy  could  be  indicative  of  serious 
problems  in  the  fluid  dynamics  model.  Dr.  P.  Crowell  and  co-workers  at  Logicon  RDA 
(NM)  speculated  that  condensation  of  H20  in  the  supersonic  expansion  might  be 
responsible  for  the  downstream  heating.  Inclusion  of  this  process  in  their  model  could 
reproduce  the  observed  temperatures,  but  only  when  it  was  assumed  that  all  of  the  H20  in 
the  flow  condensed.  The  groups  at  both  Kirtland  AFB  and  Logicon  RDA  requested 
experimental  investigation  of  this  problem. 

We  searched  for  the  condensation  heating  effect  by  observing  the  temperature  in  a 
He/l2/H20  supersonic  expansion  as  a  function  of  added  H20.  The  2-D  Laval  nozzle 
apparatus  described  in  reference  1  was  used  for  these  experiments.  Two  sample  bombs, 
each  filled  with  about  40ml  of  distilled  water,  and  a  sample  cell  containing  solid  I2  were 
added  to  the  inlet  system.  The  sample  holders  were  configured  to  pass  the  helium  over  the 
HjO  and  then  entrain  I2.  For  measurements  of  Ij/He  alone,  water  vapor  was  removed  from 
the  flow  by  either  bypassing  the  sample  bombs  or  using  a  dry  ice/acetone  mixture  to  trap 
the  water  vapor.  Both  methods  proved  equally  effective.  The  amount  of  water  vapor  added 
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was  varied  between  0  and  300  mTorr  (downstream),  as  determined  by  monitoring  the 
increase  in  the  pressure  of  the  expansion. 

The  temperature  of  the  expansion  with  and  without  the  addition  of  water  vapor  was 
estimated  spectroscopically  by  using  a  dye  laser  to  excite  the  B-X  (14-0,  17-1)  transitions  of 
I2.  The  relative  intensities  of  the  rotational  and  vibrational  transitions  were  simulated  to 
estimate  the  local  temperature  of  the  expansion.  The  temperature  was  determined  1  cm  and  4 
cm  downstream  from  the  throat  of  the  nozzle.  Typically  the  plenum  pressure  was  around  10 
Torr,  resulting  in  a  downstream  pressure  of  about  2  Torr.  Due  to  the  limitations  of  our 
pumps,  we  were  unable  to  use  plenum  pressures  much  above  10  Torr.  Using  the  highest 
B  O  partial  pressure  we  could  achieve  (300  mTorr  downstream)  and  the  maximum  plenum 
pressure,  we  were  unable  to  detect  any  heating  effects  due  to  the  presence  of  BjO. 

As  an  alternative  means  of  explaining  the  high  temperature,  it  was  proposed  that 
quenching  of  I*  by  HjO 

I*  +  Hp  — >  1+  Up  {48 }‘ 

might  be  responsible  for  the  heat  release.  However,  the  rate  constant  for  this  process  would 
have  to  be  greater  than  the  accepted  value  for  this  explanation  to  be  viable.  At  the  request  of 
Dr.  Gordon  Hager  (AFRL,  Kirtland  AFB)  and  Dr.  Peter  Crowell  we  re-examined  the  rate 
constant  for  reaction  {48 } .  Measurements  were  made  at  room  temperature  under  slow  flow 
conditions.  We  obtained  a  quenching  rate  constant  of  k6=  (2.4±0.5)xl012  cmW,  in 
excellent  agreement  with  the  results  of  several  previous  studies. 

We  attempted  to  explore  the  possibility  that  the  rate  constant  for  reaction  {48} 
might  increase  with  decreasing  temperature.  Unfortunately,  it  proved  to  be  very  difficult  to 
obtain  stable  H20  concentrations  in  the  Laval  expansion,  and  we  were  unable  to  make  a 
reliable  rate  constant  determination  at  150  K.  However,  the  results  were  not  consistent  with 
a  dramatic  increase  in  the  rate  constant  in  going  from  300  to  150  K.  It  is  unlikely  that 
reaction  {48}  can  account  for  the  relatively  high  temperatures  seen  in  VertiCOIL. 

As  the  models  seem  to  be  unable  to  account  for  the  heat  release,  the  Kirtland  group 
is  carefully  re-evaluating  their  temperature  measurements.  The  temperatures  were 
determined  from  the  lineshape  of  the  I*-I  3-4  transition5.  This  line  was  observed  using  the 
absorption  of  light  from  a  tunable  diode  laser.  Turbulence  in  the  gas  flow  near  the  walls 
and  observation  ports  could  influence  the  measurements  and  lead  to  overestimation  of  the 
temperature. 
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4.  Mechanism  of  I2  dissociation  in  COIL  systems 

The  mechanism  by  which  02(a)  dissociates  I2  in  COIL  systems  is  a  long-standing 
unsolved  problem6.  This  is  an  important  issue  as  the  dissociation  process  significantly 
impacts  the  efficiency  of  the  laser.  There  are  two  mechanisms  by  which  this  occurs.  First, 
as  it  takes  an  average  of  4  to  6  02(a)  molecules  to  dissociate  one  I2  molecule,  the 
dissociation  process  consumes  an  appreciable  amount  of  energy.  Second,  the  time  scale  of 
the  dissociation  process  influences  energy  extraction  as  it  determines  the  downstream 
position  at  which  the  maximum  concentration  of  I*  will  occur.  The  efficiency  and  time 
scale  of  the  dissociation  process  depend  on  concentration  ratios,  mixing  dynamics  and  total 
pressure.  At  present  our  ability  to  design  new  types  of  mixing  nozzles  and/or  develop 
systems  that  operate  at  higher  pressures  is  limited,  in  part,  by  our  poor  understanding  of  the 
dissociation  kinetics. 

The  currently  accepted  model  for  I2  dissociation  is  derived  from  flow  tube  studies 
conducted  by  Heidner  et  al.21  These  investigators  were  unable  to  find  a  unique  kinetic 
model  to  account  for  their  data.  They  offered  two  limiting  models  that  defined  a  range  of 
rate  constants  for  the  elementary  reaction  steps.  The  rate  constant  package  used  to  predict 
COIL  performance  is  mostly  based  on  one  of  the  limiting  models19*20.  The  mechanism  for 
Ij  dissociation  assumes  that  vibrationally  excited  L,(X)  is  the  immediate  precursor  of  atomic 
Iodine19-20.  However,  subsequent  measurements  of  I^X)  vibrational  relaxation  rate 
constants  yielded  values  that  were  not  in  agreement  with  the  rate  constants  used  in  COIL 
models6-22. 

As  there  is  appreciable  flexibility  in  the  kinetic  model,  we  considered  the  possibility 
that  improvements  could  be  made  by  constraining  the  vibrational  relaxation  rate  constants 
and  treating  some  of  the  less  well  known  rate  constants  as  variable  parameters7.  The 
original  flow  tube  data  of  Heidner  et  al.21  were  used  in  this  re-optimization  procedure.  This 
exercise  was  carried  out  in  collaboration  with  Capt.  J.  S.  Paschkewitz  (Computational 
Sciences  Branch,  Air  Vehicles  Directorate,  AFRL,  Wright  Patterson  AFB).  Improved  rate 
constants  were  obtained  that  are  significantly  different  from  those  used  in  current 
computational  models.  Unfortunately,  although  the  new  rate  constants  provided  a  better 
representation  of  the  experimental  data,  troubling  discrepancies  remained.  To  make  further 
progress  we  examined  the  possible  participation  of  electronically  excited  states  of  I2 
(denoted  by  I2(A  ’,A))8.  Channels  involving  I2(A’,A)  were  added  to  the  standard  COIL 
model  and  attempts  were  made  to  re-optimize  the  rate  constant  package  by  fitting  the  I2 
dissociation  rate  data.  In  the  initial  study  the  rate  constant  for  the  reaction 

*  In  sections  1  and  2  the  numbering  of  reactions  is  chosen  to  be  consistent  with  the  COIL  standard  rate 
constant  package19-20.  The  number  for  each  reaction  is  given  in  curly  brackets. 
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l^X,  v>20)  +  02(a‘A)  ->  21  +  02(X3Z)  {34} 

was  not  reduced  below  k^l.SxlO'11  cm3  s'1  (1/20  the  value  used  in  the  standard  package). 
With  such  rapid  dissociation  of  I^X,  v>20)  the  channel  through  the  electronically  excited 
states  could  not  make  a  significant  contribution.  Preliminary  indications  were  that  the 
deficiencies  of  the  standard  model  could  not  be  overcome  by  including  the  participation 
electronically  excited  I27.  However,  restricting  the  range  of  values  considered  for  k34  was 
somewhat  arbitrary  as  there  are  no  experimental  data  for  this  reaction.  Greater  reduction  of 
forces  the  dissociation  process  to  proceed  via  electronic  excitation  of  I2.  Under  these 
conditions  the  kinetics  become  sensitive  to  electronic  quenching  of  1 2(A',A).  A  model  of 
this  form  has  the  potential  to  give  a  good  representation  of  the  dissociation  rate  data, 
provided  that  quenching  by  Ar  is  very  slow  and  that  H20  is  more  effective  than  02(X)  in 
dissipating  energy  from  I2(A',A).  It  is  probable  that  both  the  A  and  A  states  of  I2  exhibit 
similar  quenching  kinetics.  Hence,  for  simplicity  the  remainder  of  this  discussion  is 
focused  on  I2(A0  (in  addition,  this  species  is  appreciably  more  metastable  than  I2(A)). 
Tellinghuisen  and  Phillips23  have  shown  that  quenching  of  \2{Ar)  by  Ar  is  slow  (k=2.8xl014 
cm3  s'1)  and  that  the  primary  mechanism  is  collision  induced  dissociation.  This  seemed 
encouraging,  but  data  for  quenching  by  02(X)  and  HjO  were  needed  to  further  examine  the 
role  of  I2(A')  in  the  dissociation  process. 

As  the  required  rate  constants  for  quenching  of  I^AO  by  02  and  H20  were  not 
available  we  measured  these  parameters  using  pulsed  laser  pump-probe  techniques24.  To 
populate  the  A'  state  a  mixture  of  Ij  in  Ar  was  excited  by  193  nm  pulses  from  an  ArF  laser. 
This  initially  promotes  I2  to  the  D(0*u)  state.  Collisions  with  Ar  rapidly  induce  transfer  to 
the  D{ 2g)  state,  which  then  radiates  down  to  A'3I1(2U).  This  sequence  occurs  within  the 
radiative  lifetimes  of  the  D  and  D'  states  (approximately  10  ns).  The  population  in  the  A' 
state  was  monitored  by  probe  laser  excitation  of  the  D’<—A'  transition.  Loss  of  population 
from  A'  was  followed  by  varying  the  delay  between  the  pump  and  probe  lasers.  For  low 
concentrations  of  I2  in  20  Torr  of  pure  Ar  the  A'  state  decayed  with  a  lifetime  of  about  30 
ps.  This  was  consistent  with  the  earlier  study  by  Tellinghuisen  and  Phillips23.  Metered 
quantities  of  02  or  H20  were  added  to  the  I/Ar  mixtures  to  observe  quenching  by  these 
collision  partners.  The  quenching  rate  constants  were  found  to  be  (6.3±0.6)xl0‘12  and 
(3.4±0.4)xl012  cm3  s'1  for  02(X)  and  HjO,  respectively. 

These  measurements  showed  that  both  02  and  H20  are  much  more  effective 
quenchers  of  I^A1)  than  Ar.  However,  Heidner  et  al.’s21  models  require  H20  to  be  a  better 
deactivator  than  Oz  whereas  the  present  results  indicate  the  opposite  trend.  This  apparent 
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disagreement  can  be  resolved  if  quenching  by  02  is  dominated  by  the  Franck-Condon 
favored  transfer  process 

I2(A0  +  02(X3D  ->  I2(X,  v>10)  +  02(a‘A)  { 28 } 

Data  for  matrix  isolated  ^J02  mixtures  shows  that  this  is,  indeed,  a  facile  transfer  process25. 

To  test  the  notion  that  I2(A')  is  a  significant  dissociation  intermediate,  a  new  kinetic 
model  of  the  dissociation  process  has  been  examined24.  This  model,  which  is  shown 
schematically  in  Fig.  5,  is  based  on  the  standard  reaction  set  with  the  addition  of  reactions 
involving  I2(X,  10<v<20)  andI2(A').  The  full  set  of  reactions  and  rate  constants  is  listed  in 
Table  1.  Numerical  integration  of  the  coupled  rate  equations  was  used  to  predict  I2(X) 
removal  kinetics.  Several  rate  constants  were  manually  adjusted  to  achieve  reasonable 
agreement  with  the  dissociation  rate  data  of  Heidner  et  al.21  The  revised  model  gives  a 
significantly  better  fit  to  the  experimental  data  than  the  standard  package.  From  a  chemical 
physics  perspective  the  most  appealing  aspect  of  the  model  is  that  the  rate  constants  all  have 
magnitudes  that  are  appropriate  for  the  types  of  process  that  they  are  associated  with.  Note 
that  the  specific  values  for  the  rate  constants  in  Table  1  are  not  unique  and  they  are  not  fully 
optimized.  Improvements  could  be  achieved  by  coupling  the  differential  equation  solver  with 
a  non-linear  least  squares  fitting  routine. 

The  ability  of  this  new  model  to  predict  the  performance  of  COIL  systems  is  being 
tested  by  Dr.  T.  Madden  (AFRL  Kirtland)  and  Dr.  B.  Barmashenko  (Ben  Gurion 
University,  Israel).  Preliminaiy  results  from  the  group  at  BGU  indicate  that  the  model  gives 
much  improved  predictions  for  the  dependence  of  the  laser  power  on  the  reagent 
concentrations  and  flow  conditions. 
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Table  1. 
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Figure  5.  Schematic  of  the  proposed  dissociation  mechanism 


Kinetic  scheme  used  to  simulate  the  experimental  data  of  Heidner  et  al.21 
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Reaction  #  Reaction 


Rate  Constant/  cm3  s' 


1 

02(a)+02(a)-»02(b)+02(X) 

2.7xlO'17 

21 

02(b)+I2— »02(X)+2I 

4.0x1  O'12 

40 

02(a)+I— »02(X)+I* 

7.6xl0'n 

41 

02(X)+I*— »02(a)+I 

2.6xl0'n 

44 

02(a)+I*-»02(b)+I 

l.lxlO'13 

32 

02(a)+I2— »02(X)+I2*  * 

2.0xl0'14 

34 

02(a)+I2*  *— »02(X)+2I 

lxlO'13 

27 

02(X)+I2**->02(X)+I2* 

3xl012 

39 

h2o+i2**-»h2o+i2* 

5x1  O'12 

38 

I2*+02(a)^I2(A,)+02(X) 

6xl0'12 

25 

I2(A')+02(a)->2I+02(X) 

2xlO'10 

28 

I2(A')+02(X)-^I2*+02(a) 

6.3xl0'12 

29 

I2(A')+H20^I2+H20 

3.4xl0*12 

63 

i2*+h2o->i2+h2o 

1.2xl0'12 

64 

I2*+02(X)->I2+02(X) 

1.5xl0'12 

65 

I2**+Ar-»I2*+Ar 

2.0xl0'12 

66 

I2*+Ar— >I2+Ar 

2.0xl0'13 

33 

I*+I2 — >I+I2*  * 

3.5xl0*u 

5 

02(b)+H20-»02(X)+H20 

5.5xl0'12 

48 

i*+h2o->i+h2o 

2.0xl0*12 

67 

02(b)+I2*— >02(X)+2I 

3.0xl0'12 

Boldface  type  is  used  to  indicate  reactions  that  are  not  part  of  the  standard  package  and  rate 
constants  that  have  been  assigned  values  that  are  different  to  those  used  in  the  standard 
package.  The  reaction  numbering  scheme  is  chosen  to  be  consistent  with  Table  1  of 
reference  20  l**  is  I2(X,v>20)  and  I*  is  I2(X,  10<v<20). 
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5.  Kinetic  Studies  of  the  NCI  radical. 

Jo.  Quenching  rate  constants  for  NCl(a) 

Interest  in  the  possibility  of  using  NCl(a)  as  the  energy  carrier  in  iodine  laser 
systems9'11  motivated  several  previous  studies  of  NCl(a)  removal  kinetics12-14.  Two  different 
experimental  techniques  were  used  to  determine  quenching  rate  constants,  and  the  results 
appeared  to  be  strongly  dependent  on  the  technique  used.  Hewett  et  al.14  performed 
measurements  in  a  discharge  flow  system.  NCl(a)  was  generated  from  the  reaction 

Cl  +  N3  -4  NCl(a)  +  N2  (5) 

and  the  intensity  of  the  NCI  a-X  system  was  monitored  as  a  function  of  quenching  gas 
density  and  time  under  steady  state  conditions.  Alternatively,  Ray  and  Coombe12,  and 
Henshaw  et  al.13  used  pulsed  photolysis  of  C1N3  to  generate  NCl(a),  and  time-resolved 
fluorescence  measurements  to  determine  quenching  rate  constants.  As  photolysis  of  C1N3 
produces  strong  emissions  from  the  N2  B-A  system,  an  interference  filter  was  used  to  isolate 
the  NCI  a-X  fluorescence.  Ray  and  Coombe12  reported  that  NCl(a)  generated  by  248  nm 
photolysis  of  C1N3  has  a  complex  time  profile.  Under  certain  conditions  they  found  that  the 
NCl(a)  signal  could  show  rise-fall  characteristics.  This  behavior  was  attributed  to  secondary 
generation  of  NCl(a)  from  chain  decomposition  of  C1N3.  Chain  decomposition  did  not 
appear  to  be  of  importance  when  193  nm  photolysis  was  used.  Instead,  bi-exponential 
decays  were  observed.  The  long-lived  component  of  the  decay  was  attributed  to  NCl(a). 
Ray  and  Coombe12,  and  Henshaw  et  al.13  (who  also  used  193  nm  photolysis)  determined 
NCl(a)  removal  rate  constants  by  observing  the  dependence  of  the  long-lived  decay  rate  on 
the  pressure  of  added  quenchers  or  reactants.  The  level  of  agreement  between  the  rate 
constants  reported  by  these  two  groups  was  good,  but  several  of  the  constants  were 
dramatically  different  from  values  measured  in  a  discharge  flow  system  by  Hewett  et  al.14 
We  re-examined  the  photolysis  approach  to  probe  the  origin  of  these  discrepancies  and 
determine  which  data  are  most  reliable. 

In  our  experiments  NCl(a)  was  generated  by  photolysis  of  C1N3.  The  apparatus  used 
for  these  measurements  is  shown  in  Fig.  6.  The  generator  used  to  provide  C1N3  was  the 
same  as  that  described  by  Henshaw  et  al.13  The  azide  was  formed  when  a  mixture  of  5%  Cl2 
in  He  was  passed  over  the  surface  of  moist  NaN3.  A  standard  drying  agent  was  used  to 
remove  water  from  the  CINj/Clj/He  mixture  at  the  exit  of  the  generator.  The  concentration 
of  CINj  was  monitored  via  the  absorption  band  at  250  nm.  After  the  generator  had  been  run 
for  a  few  hours,  the  efficiency  for  converting  Cl2  to  C1N3  was  typically  85-95%.  Due  to 
incomplete  conversion,  a  background  concentration  of  Cl2  was  always  present  in  the 
photolysis  cell. 
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Figure  6.  Apparatus  used  to  observe  the  fluorescence  decay  and  quenching  of  NCI  {a). 
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Figure  7.  Fluorescence  decay  curves  for  NCl(a)  in  the  presence  of  02 
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Low  photolysis  laser  powers  were  used  to  photodisociate  the  C1N3  for  the  NCl(a) 
quenching  rate  constant  measurements.  Emissions  resulting  from  photolysis  were  dispersed 
by  a  0.2  m  monochromator.  For  detection  of  fluorescence  with  X>700  nm,  a  695  nm  long- 
pass  filter  was  placed  between  the  photolysis  cell  and  the  monochromator  to  block  scattered 
laser  light  and  short  wavelength  emission  signals.  The  dispersed  light  was  detected  by  a 
cooled  photomultiplier  tube  (Hamamatsu  R1767,  SI  photocathode).  To  record  emission 
spectra,  the  output  from  the  phototube  was  processed  by  a  boxcar  integrator  (SRS  250).  The 
time  dependence  of  the  fluorescence  signals  was  examined  using  a  digital  storage 
oscilloscope  (LeCroy  ScopeStation  140, 100MHz  bandwidth).  Signal  averaging  was  used  to 
improve  the  signal  to  noise  ratios. 

For  193  nm  photolysis  of  C1N3  we  found  that  the  emission  spectrum  in  the  vicinity 
of  the  NCI  a-X  0-0  band  was  very  congested.  The  interfering  emissions  included  strong 
bands  from  the  N2  B-A  system.  Attempts  to  determine  NCl(a)  removal  rate  constants  using 
193  nm  photolysis  produced  very  scattered  results.  Surprisingly,  we  found  that  248  nm 
photolysis  produced  an  emission  spectrum  dominated  by  the  NCI  a-X  0-0  band.  The  time 
dependence  of  this  emission  gave  a  good  fit  to  a  single  exponential  decay.  Our  initial 
observations  were  made  using  the  unfocused  output  from  an  excimer  laser,  which  provided 
fluences  of  the  order  of  30  mJ  cm2  in  the  cell.  The  extent  of  chain  decomposition  should 
be  dependent  on  the  laser  power,  and  may  be  insignificant  when  low  fluences  are  used. 
However,  when  the  laser  was  focused  by  a  1  m  focal  length  lens  we  were  still  unable  to 
observe  temporal  profiles  that  showed  evidence  of  secondary  NCl(a)  production. 

As  248  nm  photolysis  appeared  to  be  suitable  for  kinetic  measurements  in  our 
apparatus,  we  examined  removal  of  NCl(a)  by  02,  Cl2,  HC1,  and  H215.  Typical  fluorescence 
decay  curves  for  various  pressures  of  02  are  shown  in  Fig.  7.  Stem-Volmer  plots  for 
removal  by  02,  Cl2,  andHCl  are  shown  in  Fig.  8.  Quenching  rate  constants  determined  in 
the  present  study  are  collected  in  Table  2,  where  they  may  be  compared  with  the  results  of 
previous  measurements.  The  rate  constant  for  quenching  by  02  was  greater  than  previous 
values  by  a  factor  of  two.  However,  the  higher  value  was  supported  by  the  time-resolved 
data  for  NC1(X)  (described  below).  The  rate  constants  for  removal  by  Cl2,  HC1,  and  H2 
were  much  closer  to  the  flow  tube  results14  than  the  values  obtained  using  193  nm 
photolysis12,13.  At  present  we  do  not  know  the  reason  for  the  large  discrepancy  between  the 
results  of  Hewett  et  al.14  and  those  of  Ray  and  Coombe12,  and  Henshaw  et  al.13 
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Figure  8.  Stem-Volmer  plots  for  quenching  of  NCl(a)  by  02,  Cl2,  and  HC1. 


Table  1 

Comparison  of  Rate  Constants  for  Removal  of  NCl(a) 


Collision  partner 

1 

Rate  Constant  /10' 

2 

14  cm3  s' 

3 

4 

o2 

560±40 

280±60 

250+20 

— 

Cl2 

59±8 

40±10 

1800±300 

2900±600 

HC1 

1.8±0.4 

1.5±0.3 

490±70 

-- 

h2 

<0.1 

<0.1 

68±7 

-- 

N3C1 

1.  Present  work 

30±10 

— 

— 

57±4 

2.  K.  B.  Hewett,  G.  C.  Manke,  D.  W.  Setser,  and  G.  Brewood,  J.  Phys.  Chem.  A, 

104, 539-551  (2000). 

3.  A.  J.  Ray  and  R.  D.  Coombe,  J.  Phys.  Chem.  98,  8940  (1994) 

4.  T.  L.  Henshaw,  S.  D.  Herrera,  G.  W.  Haggquist,  L.  A.  Schlie,  J.  Phys.  Chem.  A.  101, 
4048  (1997). 
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5b.  Removal  of  NCl(a)  by  I2. 

Quenching  of  NCl(a)  by  I2  is  described  separately,  as  the  characteristics  of  this 
process  were  somewhat  different  to  those  of  the  quenching  agents  described  above.  Stem- 
Volmer  plots  for  NCl(a)+I2  were  linear,  but  the  slopes  were  dependent  on  the  repetition 
frequency  of  the  photolysis  laser.  The  apparent  quenching  rate  constant  increased  from 
around  6x10  "  cm3  s'*  at  10  Hz  to  l.OxlO10  at  0.5  Hz.  Further  increases  in  the  rate 
constant  did  not  occur  for  repetition  rates  slower  than  0.5  Hz.  These  data  indicated  that 
NCl(a)  was  being  removed  by  reaction  with  and  that  the  gas  sample  was  not  being  fully 
replenished  between  laser  shots  for  repetition  frequencies  above  0.5  Hz.  To  avoid  this 
problem,  measurements  on  the  C1N3/I2  system  were  carried  out  using  a  pulse  repetition 
frequency  of  0.2  Hz.  Under  these  conditions  the  removal  rate  constant  for  I2  was 
(1.0±0.3)xl010  cm3  s'1.  Note  that  the  rate  constant  at  10  Hz  was  still  relatively  large,  which 
suggests  that  NCl(a)  is  efficiently  quenched  by  one  or  more  of  the  reaction  products. 

To  test  the  hypothesis  that  NCl(a)  was  being  removed  by  a  reactive  process,  we 
monitored  the  time-dependent  behavior  of  I2  following  photolysis  of  He/I/CINj  mixtures. 
Fig.  9  shows  the  apparatus  used  for  these  measurements.  The  488  nm  radiation  from  an  Ar 
ion  laser  was  used  to  monitor  I2  via  the  B-X  continuum  absorption.  An  advantage  of  this 
approach  is  that  the  strength  of  the  absorption  signal  is  insensitive  to  vibrational  excitation 
of  ^(X).  The  time  dependent  behavior  of  the  absorption  signal  is  shown  by  trace  (a)  in  Fig. 
10.  The  initial  negative-going  spike  (0-200  ps)  in  this  trace  was  caused  by  radio  frequency 
noise  from  the  pulsed  laser.  The  absorption  signal  recovered  with  a  time  constant  of  about 
340  s'1.  This  was  consistent  with  diffusion  of  I2  into  the  photolysis  region  from  the 
surrounding  un-irradiated  gas.  Trace  (b)  was  recorded  with  just  an  IfHe  mixture  present  in 
the  cell.  This  measurement,  which  shows  only  the  initial  noise  spike,  was  made  to 
demonstrate  that  ^  was  not  being  dissociated  by  the  248  nm  pulse. 

Photolysis  of  He/CIN^  mixtures  yielded  emission  from  I*.  Fig.  11a  shows  the 
time  dependence  of  the  I*  signal  for  various  initial  I2  pressures.  Note  that  I*  signals  were 
not  observed  for  248  nm  excitation  of  flowing  mixtures  of  He/I2  or  He/Clj/Ij.  The  finite 
rise  time  of  the  I*  signal  for  photolysis  of  He/C1N3/I2  mixtures  indicated  that  I*  was  being 
generated  by  secondary  photochemical  reactions.  Given  that  NCl(a)  reacts  with  I2  it  is 
probable  that  I*  was  generated  by  the  sequence 


C1N3  +  hv  NCl(a)  +  N2 

(6) 

NCl(a)  +  I2->  INC1  + 1 

(7) 

NCl(a)  + 1  ->  NC1(X)  + 1* 

(8) 
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I*  Fluorescence  Intensity 


Figure  11a.  I*  emission  following  photolysis  of  He/C1N3/I2  mixtures 


Figure  lib.  Kinetic  simulations  of  the  curves  shown  in  Fig.  11a. 
See  text  for  details. 
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We  simulated  the  I*  temporal  profiles  using  a  kinetic  model  that  included  the  above 
reactions  (6, 7,  and  8)  and  the  additional  processes 


NCl(a)  +  C1N3  -4  products 

(9a) 

NCl(a)  +  C1N3  -4  NC1(X)  +C1N3 

(9b) 

NCl(a)  +  Cl2  -4  NC12  +  Cl 

(10) 

I*  +  I2  -4  I  + 12 

{33} 

I*  +  C1N3  -4  products 

(ID 

I*  +  INC1  -4  I  +  INC1 

(12) 

NCl(a)  +  INC1  -4  NC1(X)  +  INC1 

(13) 

Literature  values  were  available  for  the  rate  constants  needed  for  this  model  with  the 
exception  of  the  constants  for  reactions  12  and  13.  However,  attempts  to  simulate  the  I* 
kinetics  (Fig  11a)  were  unsuccessful  regardless  of  the  values  chosen  for  kI2  and  k13.  The 
problem  was  traced  to  the  quenching  of  I*  by  C1N3,  which  proceeded  much  too  rapidly  in 
the  model.  To  obtain  agreement  with  the  experimental  data  the  rate  constant  for  reaction  1 1 
had  to  be  reduced  to  values  below  2xl0'14  cm3  s'1.  The  model  also  required  relatively  large 
rate  constants  for  quenching  of  I*  and  NCl(a)  by  the  product(s)  of  reaction  7.  The  best 
simulation,  which  is  shown  in  Fig.  lib,  was  obtained  using  k12=4xlO'H  and  k13=2.5xl0 11 

3  -1 

cm  s  . 

In  the  above  model  we  have  tentatively  assumed  that  the  reaction  of  NCl(a)  with  I2 
gives  INC1  as  a  product.  This  is  analogous  to  the  reaction  between  NCl(a)  and  Cl2.  To  date 
there  have  been  no  reported  observations  of  INC1,  but  ab  initio  calculations  predict  that  this 
radical  is  bound  by  at  least  1  eV. 

5c.  Detection  and  characterization  of  NCl(X)  via  b-X  absorption  spectroscopy 
(collaboration  with  Dr.  G.  C.  Manke  II  (AFRL,  NM)  and  Dr.  S.  J.  Davis(PSI,  MA)). 

A  method  for  detecting  NC1(X)  via  the  b-X  absorption  bands  was  developed  for 
three  reasons.  First,  the  b-X  bands  provide  a  convenient  diagnostic  for  the  conditions  within 
an  NCI/I  laser  system.  High-resolution  absorption  measurements  can  provide  local 
temperature  and  concentration  information.  Second,  time-resolved  monitoring  of  NC1(X)  in 
kinetics  experiments  can  be  used  to  learn  more  about  the  collision  dynamics  of  both  NC1(X) 
and  NCl(a).  Thild,  monitoring  of  NC1(X)  was  used  to  examine  the  NCl(a)/NCl(X)  ratio 
produced  by  photolysis  of  C1N3. 

The  apparatus  used  to  observe  the  NC1(X)  absorption  spectra  is  shown  in  Fig.  12. 
NCI  was  generated  by  photolysis  of  C1N3  (at  193  or  248  nm).  The  output  from  a  CW  ring 
dye  laser  was  used  to  observe  the  b-X  bands  near  665  nm. 
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Figure  13.  Section  of  the  NCI  b-X  absorption  spectrum 


29 


The  beam  was  split  into  reference  and  sample  beams.  White  cell  mirrors  were  used  to 
obtain  about  100  passes  of  the  sample  beam  through  the  absorption  cell.  The  sample  and 
reference  beams  were  then  processed  by  a  noise  canceling  detector  system  (New  Focus 
Nirvana  detector). 

Atypical  NCI  absorption  spectrum  from  this  apparatus  is  shown  in  Fig.  13.  Note 
that  the  signal-to-noise  ratio  of  this  trace  was  limited  by  the  scanning  electronics  that  control 
the  dye  laser.  The  minimum  sweep  speed  was  not  slow  enough  to  permit  extensive  signal 
averaging.  The  relative  positions  of  the  b-X  0-0  band  rotational  lines  woe  in  excellent 
agreement  with  the  work  of  Colin  and  Jones26.  However,  we  found  a  small  systematic 
calibration  error  in  the  former  study.  Our  measurements  indicate  that  the  line  positions  of 
Colin  and  Jones  should  be  reduced  by  0.04  cm'1.  This  correction  is  of  significance  for  the 
development  of  diode  laser  diagnostic  instruments.  Fig.  13  and  similar  spectra  were 
recorded  in  the  presence  of  6  Torr  of  buffer  gas  (mostly  He),  with  a  100  jxs  delay  between 
the  photolysis  pulse  and  the  absorption  measurement  (a  boxcar  integrator  was  used  to 
monitor  absorption  after  a  fixed  delay).  Under  these  conditions  the  rotational  and 
translational  energy  of  NCI  should  be  in  thermal  equilibrium  with  the  buffer  gas.  This  was 
confirmed  by  the  rotational  line  intensity  distribution,  which  defined  a  rotational  temperature 
of  300±20  K.  The  translational  temperature  was  estimated  from  absorption  line  profiles. 
The  lineshapes  were  well  represented  by  a  Gaussian  function  with  a  linewidth  of  0.031  cm'1 
(FWHM).  This  is  slightly  wider  than  would  be  expected  for  a  temperature  of  300  K 
(predicted  width  of  0.026  cm'1).  However,  fitting  to  a  single  Gaussian  function  ignores  the 
hyperfine  structure  of  the  ro-vibronic  transitions.  We  have  calculated  the  hyperfine 
structure  of  the  line,  and  found  that  the  transitions  span  a  range  of  0.05  cm'1.  Convolution 
of  the  hyperfine  structure  with  the  Doppler  lineshape  for  300  K  readily  accounts  for  the 
observed  linewidth.  Clearly,  diagnostic  techniques  that  use  the  profiles  of  NCI  rotational 
lines  for  temperature  determination  will  need  to  take  the  hyperfine  structure  of  the  fines  into 
account. 

Searches  were  made  for  transitions  originating  from  vibrationally  excited  levels  of 
NC1(X),  to  examine  the  possibility  that  vibrationally  excited  products  are  produced  by 
photolysis  of  C1N3.  We  were  unable  to  detect  population  in  the  v=l  or  2  levels  using  either 
193  or  248  nm  photolysis.  From  these  measurements,  and  estimates  of  the  sensitivity  of  the 
apparatus,  we  conclude  that  60%  or  more  of  the  NC1(X)  is  produced  in  the  v=0  level. 

The  primary  product  resulting  from  UV  photolysis  of  C1N3  is  NCl(a).  Previous 
studies27  have  shown  that  02  quenches  NCl(a)  by  the  E-E  transfer  process 

NCl(a)  +  02(X)  -4  NC1(X)  +  02(a)  (14) 
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In  the  absorption  measurements,  addition  of  02  to  the  photolysis  cell  produced  an  increase 
in  the  NC1(X)  concentration.  For  example,  Fig.  14  shows  scans  over  the  ^(14)  line  with 
and  without  02  present  in  the  photolysis  cell.  These  traces  were  recorded  using  193  nm 
photolysis  light.  The  dependence  of  the  NC1(X)  signal  on  02  pressure  was  consistent  with 
a  competition  between  quenching  and  removal  of  NCl(a)  by  processes  that  do  not  generate 
NC1(X)  (e.g.,  NCl(a)+ClN3  -4  NC12  +N3).  For  02  partial  pressures  above  0.5  Torr  the 
NC1(X)  signal  was  not  perceptibly  dependent  on  further  addition  of  02,  indicating  that 
NCl(a)  removal  was  dominated  by  02  quenching  under  these  conditions.  If  it  is  assumed 
that  all  of  the  NCl(a)  created  by  photolysis  is  quenched  by  02  at  the  high  pressure  limit, 
data  such  as  those  shown  in  Fig.  14  yield  a  branching  fraction  for  N3C1  +  193  nm  -4 
NCl(a)  +  N2  of  approximately  0.75.  Similar  measurements  performed  using  248  nm 
photolysis  indicate  that  the  branching  fraction  for  N3C1  +  248  nm  -4  NCl(a)  +  N2  is 
greater  than  0.8. 

Fig.  15  shows  a  typical  time-resolved  absorption  signal.  Such  data  were  recorded 
for  a  variety  of  total  pressures  and  added  quenching  gas  conditions.  The  results  were 
analyzed  by  kinetic  modeling.  For  example,  the  data  for  energy  transfer  induced  by  02 


were  analyzed  by  a  model  based  on  the  reactions 

C1N3  +  248  nm  -4  NCl(a)  +  N2  (6) 

NCl(a)  +  C1N3  -4  products  (9) 

NCl(a)  +  Cl2  -4  NC12  +  Cl  (10) 

NCl(a)  +  02(X)  -4  NC1(X)  +  02(d)  (14) 

NC1(X)  +  NC1(X)  -4  N2  +  2C1  (15) 

NC1(X)  +  C1N3  -4  products  (16) 


The  rate  constants  for  most  of  these  reactions  were  determined  from  quenching  studies  (see 
Table  2).  Clyne  et  al.28  reported  a  rate  constant  of  8.1xl0'12  cm3  s'1  for  reaction  15,  and  the 
pseudo  first  order  removal  of  NC1(X)  by  C1N3  was  represented  by  a  decay  rate  of  2800  s'1. 
The  smooth  curve  in  Fig.  15  illustrates  the  ability  of  this  kinetic  model  to  fit  the  absorption 
data.  Note  that  we  have  assumed  a  branching  fraction  for  NCl(a)  formation  of  one  for  this 
model.  In  earlier  analyses  of  these  data  we  did  not  include  the  NC1(X)+NC1(X)  reaction15. 
Without  this  process  the  kinetic  model  was  unable  to  reproduce  the  relatively  fast  rise  times 
of  the  absorption  signals.  To  account  for  the  fast  rise  times  we  assumed  that  vibrationally 
excited  NC1(X,  v>0)  was  formed  and  rapidly  relaxed  by  collisions  with  the  buffer  gas.  With 
reaction  15  in  the  model  this  assumption  was  not  required. 
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Figure  14.  Scans  over  the  NC1(X)  ^(14)  line  for  photolysis  of  He/C1N3  and  He/ClN3/02 
mixtures. 


Figure  15.  Time  resolved  absorption  of  NC1(X)  from  photolysis  of  C1N3  in  the  presence  of 
02  ([O2]=2.2x1015  cm'3).  The  smooth  curve  shows  a  fit  of  the  kinetic  model  to  the  data. 
See  text  for  details.. 
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5d.  Determination  of  the  NCl(a)  self-annihilation  rate  constant 

Self-annihilation  of  NCl(a), 

NCl(a)  +  NCl(a)  — >  products  (17) 

is  one  of  the  most  critical  reactions  for  any  laser  that  relies  on  this  metastable  as  an  energy 
carrier.  The  room  temperature  rate  constant  for  NCl(a)  self-annihilation  was  investigated  by 
Henshaw  et  al.13  who  reported  a  value  of  (7.2±0.9)xl0'12  cm3  s*1.  However,  our  work  on 
NCl(a)  decay  kinetics  indicated  that  there  may  have  been  significant  problems  with  this 
measurement15.  Consequently,  we  re-examined  NCl(a)  self-annihilation. 

The  self-annihilation  process  cannot  be  observed  under  pseudo  first-order 
conditions.  To  extract  the  rate  constant  from  second  order  decay  data,  the  absolute 
concentration  of  NCl(a)  must  be  known.  In  the  experiments  described  here  we  used  248 
nm  photolysis  of  C1N3  to  generate  NCl(a).  The  absolute  concentration  was  derived  using 
the  branching  ratio  for  formation  of  NCl(a)  established  by  the  measurements  described  in 
the  previous  section. 

Our  studies  of  NCl(a)  quenching  kinetics  were  carried  out  using  low  photolysis 
laser  fluences.  Small  concentrations  of  NCl(a)  were  formed,  and  the  fluorescence  decay 
curves  showed  no  evidence  of  second  order  removal  processes.  The  conditions  of  these 
experiments  were  such  that  initial  NCl(a)  densities  of  about  4x1 014  cm'3  were  generated.  To 
achieve  conditions  where  second  order  removal  of  NCl(a)  could  be  observed  we  focused  the 
beam  from  the  photolysis  laser  to  dissociate  a  larger  fraction  of  the  C1N3.  A  lm  focal 
length  lens  was  used  to  focus  the  beam,  and  neutral  density  filters  were  used  to  vary  the 
power  transmitted  to  the  photolysis  cell.  Imaging  optics  were  used  to  collect  the 
fluorescence  emitted  from  the  focal  volume.  The  fluorescence  was  dispersed  by  a  0.25  m 
monochromator  and  detected  by  a  photomultiplier  with  an  extended  SI  photocathode. 

Fig.’s  16  and  17  show  NCl(a)  fluorescence  decay  curves  recorded  for  a  range  of 
laser  powers.  These  traces  were  obtained  with  7  Torr  of  the  CIN/Cl/He  mixture  present  in 
the  photolysis  cell.  The  C1N3  absorption  diagnostic  indicated  that  >90%  of  the  Cl2  had  been 
converted  to  C1N3.  The  curves  in  Fig.  2  show  an  increase  in  the  initial  decay  rate  with 
increasing  laser  intensity.  This  behavior  is  readily  understood  using  the  kinetic  model 
proposed  by  Henshaw  et  al.13  Removal  of  NCl(a)  can  be  considered  in  terms  of  the  second 
order  self-annihilation  and  an  effective  rate  for  all  first-order  loss  processes  (kL),  such  that 

_d[NCl(g)]  =  ki[NC1(^j  2  +  kL[NCi(<7)]  (18) 

dt 

where  k^k^ClNJ+knfCy+kp,  and  kD  is  the  diffusion  rate.  From  this  expression  it  is 
apparent  that  the  initial  decay  rate  will  depend  on  [NCl(a)]0,  which  in  turn  depends  on  the 
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laser  intensity.  Fig.  17  shows  that,  above  a  certain  threshold,  the  NCl(a)  decay  curves 
became  independent  of  laser  intensity.  We  attribute  this  behavior  to  complete  photolytic 
destruction  of  C1N3.  Once  this  threshold  has  been  crossed  the  initial  NCl(a)  concentration 
is  no  longer  dependent  on  the  laser  intensity.  To  check  the  validity  of  this  interpretation  we 
calculated  the  intensity  needed  to  achieve  complete  dissociation  from  the  C1N3  absorption 
cross  section13  (1.5xl0*18  cm2  at  248  nm)  and  compared  this  with  the  observed  threshold 
intensity.  The  theoretical  threshold  for  complete  dissociation  (530  mJ  cm*2)  was  a  little 
lower  than  the  observed  value  (600  mJ  cm*2),  but  the  agreement  was  reasonably  good  given 
that  the  window  losses  were  not  accurately  characterized. 

To  determine  the  self-annihilation  rate  constant  the  decay  curves  recorded  using 
high  laser  intensities  (complete  dissociation)  were  fit  using  the  integrated  form  of  Eq.  18, 

kjNCKo)], _ 


[NCI(a)]  = 


exp(kLt)(k17[NCl(a)]„  +  kL)  -  k„[NCl(a)]0 


(19) 


The  initial  NCl(a)  concentration  was  calculated  by  assuming  90%  conversion  of  Cl2  to  C1N3, 
and  a  branching  fraction  of  0.8  for  generation  of  NCl(a)  from  C1N3  photolysis.  Fig.  18 
shows  an  example  of  the  level  of  agreement  obtained  between  the  data  and  the  fitted  curves. 
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From  analyses  of  seven  decay  curves  we  obtained  a  value  for  the  self-annihilation  rate 
constant  of  (6±l)xl013  cm3  s'1.  We  extended  this  analysis  to  include  data  taken  under 
conditions  where  C1N3  was  not  completely  dissociated.  For  this  situation  we  used 
measurements  of  the  laser  power  and  the  published  absorption  cross  section  to  predict  the 
initial  NCl(a)  concentration.  The  above  calculations  suggest  that  this  procedure  may 
overestimate  the  concentration  by  about  15%.  Analysis  of  all  of  the  curves  that  exhibited 
significant  second-order  decay  components  defined  a  value  for  k,7  of  (7±2)xl0 13  cm  s  . 
The  error  range  given  here  incorporates  the  uncertainty  in  the  branching  fraction. 

Computational  models  of  the  NCl(a)/I*  system  are  being  developed  by  Tim 
Madden  (AFRL,  Kirtland)  and  Peter  Crowel  (Logicon  RDA,  NM).  These  models  have 
been  optimized  to  reproduce  the  performance  of  the  azide  laser  flow  facility  at  Kirtland. 
Madden  has  examined  the  effects  of  changing  the  NCl(a)  self-annihilation  rate  constant  in 
his  CFD  model.  As  expected,  attempts  to  model  recent  I*  gain  measurements  with  a  low 
value  for  k,  overestimated  the  I*  density.  This  problem  could  be  corrected  by  reducing  the 
branching  fraction  for  the  energy  transfer  step  NCl(a)+I— >NC1(X)+I*,  which  is  not  a  well- 
defined  parameter  at  present. 

6.  Quenching  of  I*  by  C1N3  and  HN3. 

A  previous  study29  of  the  quenching  of  I*  by  C1N3  provided  a  rate  constant  of  2x10' 
11  cm3  s'1.  This  value  was  found  to  be  incompatible  with  the  data  obtained  by  photolyzing 
He/ClNj/I2  mixtures,  as  described  above.  To  further  investigate  this  conflict  we  measured 
the  quenching  rate  constant.  In  these  experiments  I*  was  generated  in  the  presence  of  C1N3 
by  photolysis  of  I2.  A  pulsed  dye  laser  operating  at  490  nm  was  used  to  dissociate  I2  via  the 
B-X  continuum.  Potential  complications  from  secondary  photochemical  reactions  were 
avoided  by  using  low  concentrations  of  I2  and  I*.  The  partial  pressure  of  I2  in  the 
photolysis  cell  was  held  at  15  mTorr  and  the  laser  intensity  (1.3  mJ  cm2)  was  sufficient  to 
produce  an  I*  density  of  approximately  2xl012  cm  3. 

Fluorescence  from  I*  was  selected  by  a  0.25  m  monochromator  and  detected  by  a 
cooled  Ge  device  (ADC  model  403).  This  detector  was  not  fast  enough  to  resolve  the  time 
dependence  of  the  I*  signal.  Hence,  the  time-integrated  signal  was  monitored  as  a  function 
of  added  quenching  agent.  To  obtain  quantitative  data  using  this  technique  it  is  necessary  to 
know  the  effective  I*  decay  rate  (mostly  determined  by  the  quenching  due  to  undissociated 
I2)  or  to  calibrate  the  measurement  using  a  species  that  has  a  well  known  quenching  rate 
constant.  We  used  the  latter  approach  with  02  as  the  calibration  species  (k=2.6xl0'ncm3  s'1 
at  298  K).  The  C1N3  quenching  rate  constant  was  derived  from  two  series  of  measurements 
made  with  the  I2  pressure  and  laser  power  held  constant.  First  the  I*  intensity  was 
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measured  as  a  function  of  added  C1N3.  The  measurements  were  then  repeated  with  C1N3 
replaced  by  02.  Fig.  19  shows  typical  Stem-Volmer  plots  (the  reciprocal  of  the 
fluorescence  intensity  vs.  quenching  gas  number  density)  obtained  by  this  procedure. 
Quenching  of  I*  by  C1N3  was  found  to  be  virtually  undetectable  in  these  experiments.  The 
slight  positive  slope  of  the  plot  (Fig.  19)  is  consistent  with  quenching  due  to  residual  Cl2 
from  the  C1N3  generator.  Comparison  with  the  quenching  by  02  placed  an  upper  limit  on 
the  rate  constant  for  quenching  by  C1N3  of  k16<2xl0'13.  This  upper  bound  is  considerably 
greater  than  the  rate  constant  used  to  model  the  C1N3/I2  photolysis  experiments. 
Unfortunately  we  could  not  refine  this  limit,  as  we  were  unable  to  remove  all  of  the  residual 
Cl2  from  the  C1N3  generator  or  accurately  quantify  the  amount  of  Cl2  present.  It  is  evident, 
however,  that  our  results  do  not  agree  with  the  previous  measurement.  The  origin  of  this 
large  discrepancy  is  unclear.  Ray  and  Coombe29  used  248  nm  photolysis  of  CF3I  to 
generate  I*  in  the  presence  of  C1N3.  The  quenching  rate  constant  was  obtained  from  a  plot 
of  the  I*  decay  rate  vs.  C1N3  pressure.  There  is  no  obvious  reason  why  this  approach 
should  yield  an  erroneous  value  for  the  rate  constant. 

Quenching  of  I*  by  HN3  is  of  interest  as  it  is  thought  that  residual  HN3  is  present 
in  the  downstream  region  of  the  azide  laser  flow  facility.  We  used  the  techniques  described 
above  to  determine  a  room  temperature  rate  constant  for  quenching  by  HN3  of  5.8xl0‘13 
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7.  Ab  initio  calculations  for  I+NC1  energy  transfer  (collaboration  with  Prof.  K. 
Morokuma) 

The  energy  transfer  process  NCl(a)+I  ->  NC1(X)  +  I*  is  the  pump  reaction  in  the 
NCl(a)/I  laser.  The  I*/I  branching  ratio  and  mechanism  of  transfer  are  both  issues  of 
relevance  to  azide-driven  laser  systems.  Following  our  recent  theoretical  study  of  the  related 
process3,  l+02(a)  ->  I*  +  02(X),  we  have  calculated  potential  energy  surfaces  for  the 
ground  and  low-lying  excited  states  of  I-NC1.  High-level  methods  were  used;  CASSCF 
and  CASPT2,  with  triple-zeta  polarized  basis  sets.  The  basis  included  diffuse  orbitals  and  a 
relativistic  effective  core  potential  for  iodine.  Calculations  that  included  spin-orbit  coupling 
were  carried  out  with  the  NCI  bond  distance  fixed  at  the  NCl(a)  equilibrium  separation  (3.0 
bohr).  The  results  were  very  much  like  those  obtained  for  I+02.  With  the  NCI  bond 
frozen,  there  are  no  deeply  bound  states  that  correlate  with  NCl(a)+I(2P3/2)  or 
NCl(X)+I(2P1/2).  Figure  20  shows  one-dimensional  cuts  through  the  surfaces  that  correlate 
with  these  limits.  The  angle  chosen  for  these  cuts  (130°)  gives  the  deepest  minima  and 
lowest  energy  curve  crossings.  As  in  the  case  of  I+02,  we  found  curve  crossings  below  the 
dissociation  asymptote  for  NCl(a)+I(2P3/2).  Other  curve  crossings  were  evident  at  higher 
energies.  Recent  measurements  of  the  temperature  dependence  of  the  NCl(a)+I  transfer 
rate  constant30  yielded  an  effective  activation  energy  of  about  750  cm1.  This  suggests  that 
the  higher  energy  surface  crossings  are  important  energy  transfer  channels. 

The  energy  defect  for  NCl(a)  +  I  transfer  is  relatively  large  (1677  cm1),  and  near¬ 
resonant  transfer  would  be  expected  to  produce  vibrationally  excited  NC1(X,  v=2).  If  the 
interaction  between  NCI  and  I  is  primarily  physical  (i.e.,  processes  occurring  on  potential 
energy  surfaces  that  are  not  significantly  bound,  as  suggested  by  Figure  20),  it  is  unlikely 
that  transfer  will  occur  with  simultaneous  vibrational  excitation  of  NC1(X).  For  this 
situation  the  excess  energy  would  be  released  to  translational  recoil,  which  is  a  dynamically 
unfavorable  process.  Conversely,  if  transfer  involves  the  transient  formation  of  a  bound 
INC1  complex,  the  excess  energy  may  be  imparted  to  NC1(X)  vibration  and/or  rotation. 

Using  spin-free  methods  we  have  searched  for  bound  states  that  correlate  with 
NCl(a)  + 1  by  varying  all  three  internal  coordinates.  Van  der  Waals  minima  were  located, 
but  none  of  the  potential  surfaces  supported  chemically  bound  states.  This  indicates  that 
transfer  is  mediated  by  physical  interactions.  However,  to  be  certain  of  this  result,  the 
calculations  need  to  be  extended  to  include  spin-orbit  coupling.  The  ground  state  surface 
(correlating  with  NC1(X)+I)  does  exhibit  a  deep  minimum  when  the  NCI  bond  is  allowed  to 
extend.  If  spin-orbit  coupling  is  strong  enough,  some  of  this  bonding  character  could  be 
mixed  in  to  the  excited  states.  Calculations  where  all  three  coordinates  are  varied,  with 
inclusion  of  spin-orbit  coupling,  will  be  carried  out  in  the  near  future. 
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Figure  20.  Potential  energy  curves  for  INC1 


8.  Ab  initio  calculations  for  Cl  +  N3  -»  NCI  +  N2  (collaboration  with  Prof.  K. 
Morokuma) 

The  reaction  Cl  +  N3  NCI  +  N2  is  known  to  give  substantial  yields  of  NCl(a). 
The  potential  use  of  this  reaction  in  a  chemically  driven  iodine  laser  system  is  currently 
being  evaluated.  The  branching  ratio  for  Cl  +  N3  ->  N2  +  NCl(a,  b)  (singlet  product)  vs.  Cl 
+  N3  N2  +  NC1(X)  (triplet  product)  is  a  critical  issue  for  laser  pumping  applications.  In 
recent  kinetic  studies  it  has  been  determined  that  NCl(a)  is  formed  with  a  yield  in  excess  of 
60%31.  However,  it  is  a  difficult  matter  to  make  a  direct  experimental  determination  of  the 
branching  ratio.  To  address  this  issue  we  are  performing  ab  initio  calculations  to  examine 
the  singlet/triplet  branching  ratio  from  a  theoretical  perspective.  Qualitative  features  of  the 
potential  energy  surfaces  have  been  investigated  using  the  B3LYP/6-31+G(d)  level  of 
theory.  Improved  energies  for  the  stationary  points  were  then  computed  using  CASSCF  (12 
electrons,  10  orbitals)  followed  by  CASPT2  (8  electrons,  6  orbitals).  The  basis  set  used 
was  of  valence  triple  zeta  quality  (spdf,  sp+ ).  Figure  21  shows  a  summary  of  the  results 
(with  the  singlet  and  triplet  states  denoted  by  S0  and  T,).  The  left  and  right  hand  sides  of 
this  figure  correspond  to  reactants  and  products,  respectively.  The  minimum  energy 
structure  for  N3C1  was  found  about  ~40  kcal/mol  below  the  entrance  channel.  There  is  no 
barrier  between  this  minimum  and  the  Cl  +  N3  reactants  side,  but  a  barrier  exists  between 
N3C1  and  the  products  NCl(a)+N2on  the  siglet  surface.  As  can  be  seen  in  Figure  21, 
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this  barrier,  which  corresponds  to  a  trans  structure,  lies  11  kcal/mol  below  the  Cl+N3 
entrance  channel,  but  33  kcal/mol  above  the  singlet  products.  The  triplet  surface  has  a 
shallower  minimum  than  S0,  lying  about  ~19  kcal/mol  below  Cl  +  N3.  There  is  no  bamer 
on  the  triplet  surface  between  Cl+N3  and  the  minimum  energy  structure  for  triplet  N3C1,  but 
a  small  barrier  (2  kcal/mol)  exists  between  N3C1(Tj)  and  the  products  NC1(X3X )  +  N2,  The 
minimum  point  on  the  S/T,  surface  of  crossing  (MSX)  occurs  before  the  S0  exit  barrier  but 
just  after  the  T,  exit  barrier,  and  is  located  ~1  kcal/mol  below  the  singlet  barrier.  As  the 
barriers  on  both  the  singlet  and  triplet  surfaces  lie  below  the  Cl+N3  entrance  channel,  the 
reason  for  the  preference  for  singlet  products  is  not  obvious.  Preliminary  4-dimensional 
scans  of  the  surfaces  were  performed  at  a  lower  level  of  theory  (CASSCF/cc-vtz  )  in  order 
to  investigate  their  topological  features.  Inspection  of  these  surfaces  suggests  that  the 
entrance  channel  for  reaction  on  the  singlet  surface  includes  a  greater  range  of  phase  space 
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than  the  entrance  for  the  triplet  channel.  This  difference  could  account  for  the  observed 
branching  ratios,  but  dynamical  calculations  are  needed  to  test  this  notion. 
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